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1. Introduction 



In the present paper we consider the initial value problem 



(dF \ ( dF \ 
x / x \ 2/ / y 



(1.1) 



u>(:c, y, 0) = <p(x, y), tu t (x, y, 0) = ^(x, j/), 



(1.2) 



where fll.ljl models anti-plane shear motions in nonlinear nonlocal elasticity, in terms 
of non-dimensional quantities. In fll.ip - fll.2p . w = w(x,y,t) represents the out-of- 
plane displacement, the strain energy density function F is a nonlinear function of 
|Vw| 2 = {w 2 x + Wy) for isotropic materials with F(0) = 0, and the subscripts denote 
partial derivatives. The terms with j3 in fll.ip incorporate the nonlocal effects where 



denotes convolution of /3 and u. The kernel /3(x,y) is assumed to be an integrable 



where C is a positive constant and r > 2. The aim of this paper is to establish the 
well-posedness of the initial value problem fll.ip - fll.2p . as well as the global existence 
and blow-up of solutions for a wide class of the kernel functions {3(x, y). The number r 
in fll.3p is closely related to the smoothness of (3 and, consequently, as the decay rate r 
gets larger the regularizing effect of the nonlocal behavior increases. This situation is 
clearly observed through a comparison of Theorem 13.71 and Theorem 13.81 

Although the model requires the nonlinearity to be of the isotropic form F(w 2 +Wy), 
our results also apply to the more general F(w x , w y )-type nonlinearities corresponding 
to the anisotropic case. The bulk of our work deals with the isotropic form of F but 
in a separate section we present all the necessary modifications corresponding to the 
anisotropic form. Similarly, in the model the kernel (3 is a function of the modulus 
\(x, y)\, but we do not require this restriction on f3 in our work. Basically, the approach 
presented here extends the techniques used for the one-dimensional nonlinear nonlocal 
Boussinesq-type wave equations in the previous studies [TJ |2j |3] to the two-dimensional 
wave equation given by f ll.ip . It is worthwhile observing that when is taken as 
the Dirac measure in f ll.ip . one recovers the quasilinear wave equation for anti-plane 
shear motions of the local (classical) theory of elasticity. A natural question is what 
happens if f ll.ip -( |L2l) is considered on a bounded domain. This question requires a 
careful interpretation of the convolution integral and of possible boundary conditions. 
We refer to the recent manuscript [1] where such interpretations are studied for nonlocal 
diffusion problems on bounded domains. 

The plan of the paper is as follows: In Section 2 we give a brief formulation of 
the anti-plane shearing problem of nonlocal elasticity. In Section 3 we present a local 
existence theory for solutions of the Cauchy problem f ll.ip - ( TP]) for given initial data in 
suitable Sobolev spaces. In Section 4 we prove global existence of solutions of (ll.ip - (ll.2p 




function whose Fourier transform, £2)) satisfies 



< AO < C(l + \i\ 2 r'\ for all i = (£i,6), 



(1.3) 
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assuming some positivity condition on the nonlinear function F together with enough 
smoothness on the initial data. In Section 5 we discuss finite time blow-up of solutions. 
Finally, in Section 6 we show that, with certain modifications, all the results of Sections 
2-5 are also valid for more general F(w x ,w y )-type nonlinearities. 

Throughout the paper, = = J R2 e~ tz ^u(z)dz and J r ~ 1 (u)(z) = 

(2^p J R 2 e 42 '^u(£)d£ denote the Fourier transform and inverse Fourier transform, 
respectively, where z = (x,y), £ = (^1,^2), dz = dxdy, and <i£ = d£,id£ 2 - Furthermore, 
H S (M. 2 ) denotes the L 2 Sobolev space on IR 2 . For the H s norm we use the Fourier 
transform representation = L 2 (l + |£| 2 ) s |^(£)| 2 d£ • Also, ||u|| ioo and 

indicate the L°° and L 2 norms, respectively, and (u, v) refers to the inner product of u 
and v in L 2 (R 2 ). 

2. The Model 

In this section we discuss how equation ( 11. ip can be derived to describe the propagation 
of a finite amplitude transverse wave in a nonlocally elastic medium. Before stating 
our derivation, we need to introduce the concept of nonlocal elasticity. One of the 
major drawbacks in the local theory of elasticity is that it does not include any intrinsic 
length scale and consequently does not take into account the long range forces that 
become increasingly important at small scales. As a result, the local theory of elasticity 
is incapable of predicting, for instance, (i) the dispersive nature of harmonic waves in 
crystal lattices and (ii) the boundedness of the stress field near the tip of a crack. In order 
to overcome such deficiencies various generalizations of the local theory of elasticity have 
been proposed. One such generalization is the theory of nonlocal elasticity which has 
been developed by Kroner [5], Eringen and Edelen [6], Kunin [7], Rogula [8], Eringen 
[91IT0] over the last several decades (For more recent studies on the subject of generalized 
theories of elasticity, see, for instance, [TH [12], [131 El EES EE] and references therein). 
What distinguishes the theory of nonlocal elasticity from the local theory of elasticity is 
that the stress at a point depends on the strain field at every point in the body. Although 
there has been a considerable amount of research done on small scale effects within the 
context of the theory of nonlocal elasticity, they are mostly restricted to linear models. 
Recently, in [TJ [21 [3] various Cauchy problems based on a one-dimensional nonlinear 
model of nonlocal elasticity have been studied. Here, we show how the approach in those 
studies is extended to the dynamic anti-plane shearing problem of nonlinear nonlocal 
elasticity. 

Consider an isotropic homogeneous nonlocally elastic medium. Identify a material 
point X of the medium by its rectangular Cartesian coordinates in a reference 
configuration: X = (X\, X2, X3). We assume that the reference configuration is 
unstressed. Let x(X, t) = (xi(X, t), x 2 (X, t), x 3 (X, t)) denote the position of the same 
point at time t. Then the displacement and the deformation gradient are given by 
u(X, t) = x(X,t) — X and A(X,t) = Grad x(X,t), respectively. We suppose that a 
(local) strain energy density function F(A) per unit volume of the undeformed reference 
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configuration exists, i.e., the material is (locally) hyperelastic, and that it sustains 
a nontrivial anti-plane shear motion. In the local theory of elasticity, a constitutive 
equation of the form cr = <x(A) = dF(A) /dA holds for a hyperelastic material (see 
equation (4.3.7) of [17J), where cr is the nominal stress tensor (note that some authors 
use its transpose referred to as the first Piola-Kirchhoff stress tensor). In the theory of 
nonlocal elasticity the (nonlocal) stress tensor S is related to the (local) stress tensor 
cr through the constitutive relation S = S(X, t) = J /3(|X — Y|)cr(A(Y, t))dY where 
/3(|X — Y|) is a kernel function that weights the contribution of the local stresses to the 
nonlocal stresses. In the absence of body forces, the (Lagrangian) equation of motion 
(see equation (3.4.4) of [17J) is given by p x = Div S where po is the mass density of the 
medium and a superposed dot indicates the material time derivative. The only difference 
between the equations of the local theory of elasticity and those of the nonlocal model 
presented here is due to the constitutive equations. 

Now we consider an anti-plane shear motion of the form 

Xl = X h x 2 = X 2} x 3 = X 3 + w(X 1 ,X 2 ,t) (2.1) 

for a nonlocally elastic material, where the out-of-plane displacement w is the only non- 
zero component of displacement, i.e. U\ = u 2 = and u 3 = w(X 1 , X 2 , t). We henceforth 
replace the arguments X\ and X 2 of the displacement w with x and y, respectively, and 
denote partial differentiations with subscript letters. For isotropic materials the strain 
energy density function F is a function of the three fundamental scalar invariants of 
the left Cauchy-Green matrix and for the anti-plane shear motion (12. ip it turns out 
to be a function of w 2 + w 2 alone: F = F{w 2 + w 2 ) (see, for instance, Section 4 of 
[18J). Furthermore, the equation of motion reduces to the scalar partial differential 
equation poW tt — (/? * cr i3) :r . + (P * °~23) y where o"i 3 and a 23 are the (local) shear stresses 
arising due to the anti-plane shear motion and they are given by o~\ 3 = dF/dw x and 
(T23 = dF/dwy. The computations are identical to those in the conventional formulation 
of nonlinear elasticity, provided we replace the nonlocal stress tensor with the local 
stress of conventional theory of elasticity [18J. The nonlocal behavior is represented by 
the convolution integral. Thus, without loss of generality, if we make a suitable non- 
dimensionalization of the equation of motion (see [2] for the non-dimensionalization in 
the one-dimensional case) and use the same symbols to avoid a proliferation of notation, 
or simply take the mass density to be 1, we get ( 11. ip . Equation ( 11. ip is consistent with 
that of the conventional formulation of nonlinear elasticity. In other words, when 
is taken as the Dirac measure to eliminate the nonlocal effect, ( 11. Ill reduces to the 
quasilinear wave equation governing anti-plane shear motions in the local theory of 
nonlinear elasticity (see for instance equation (7.10) of [18] or equation (2.2) of [19J). A 
list of the most commonly used one-dimensional kernel functions that satisfy the one- 
dimensional version of the condition given in (11. 3p is presented in [2]. We now present 
three examples of two-dimensional kernel functions used in the literature. 

(i) The Gaussian kernel [20]: &{x,y) = {2Tr)- 1 e~^ + ^^ 2 . We have = 
e -(£?+£f)/2_ This is a highly regularizing kernel as can be observed by the fact 
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that we can take any r in (11.31) . 



(ii) The modified Bess el function kernel [20]: (3(x,y) = (27r) _1 Kq(^/ x 2 + y 2 ) where 
Kq is the modified Bessel function of the second kind of order zero. Since 
= (! + £? + ^i) -1 ' f° r this special case we have r = 2 in (11.31) . Note 
that (3 is the Green's function for the operator (1 — A) where A denotes the two- 
dimensional Laplacian. In this case (II -ip can equivalently be written as 

, dF\ ( dF 
Wu - A% = - — + 



dw T , J _ \ dw y y y 



Letting F(s) — |s + G(s) we obtain the more familiar form 

, dG\ ( dG 

w tt - Aw - Aw tt = - — + 



ii 



(Hi) The bi-Helmholtz type kernel |21j : 



Z7T{C 1 C 2 ) 

where c\ and c 2 are real and positive constants. Since /3(£i, £2) = [1 + 7i(£i + £2) + 
72(^1 + ^I) 2 ] -1 where ji — c\ + c\ and 72 = c\c\ we have r = 4. As above, (3 is 
Green's function for the operator (1 — 71 A + 72 A 2 ). Then (II. ip becomes 

a 2 fdG\ f dG 
w tt - Aw - "iiAwtt + 72A w a = t: — + 



\dw x J x \dw y/y 

In the remainder of this paper we discuss the question of well-posedness of the Cauchy 
problem (Tl~Tj) - (lL2]) . 

3. Local Existence and Uniqueness of Solutions 

In the present section, we prove existence and uniqueness of solutions over a small time 
interval. Local well-posedness is established by converting the initial value problem 
( ll.ip -( Ti~2"]) into a system of Banach space-valued ordinary differential equations. Thus 
( ll.ip -l TOl) is formally equivalent to the system 

w t = v, w(0) = ip, (3.1) 
v t = Kw, v(0) = V> (3.2) 
where the operator K is defined as 

M'-S.+O"©,- (3 - 3) 

The Banach space X s will be defined as 

X s = {w G H S (R 2 ); w x ,w y G L°°(R 2 )}, 
endowed with the norm 

\\w\\x* = \\w\\ s + ||wx-||l°° + (3-4) 
The following two lemmas are useful in the proof. 
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Lemma 3.1 (Sobolev Embedding Theorem) If s > — + k, then 
H s (R n ) C C k (R n ) n L°°(R n ). 

In particular when n = 2 and s > 2, the embedding in Lemma 13.11 implies the norm 
estimate || \Vu\ < C||w|| s . We refer to Chapter 5 of [22] for a discussion on the 
many versions of the Sobolev embedding theorem. 

Lemma 3.2 Let s > and let Ui,u 2 G H s (R n ) n T/jen u x u 2 e H s (R n ) and 

||«i« 2 ||s < C , (||ui|| s ||w 2 ||l°° + ||ni|| L oo||n 2 || s ). 

Lemma [3.21 can be found in [22] in a more general L p -setting (see Lemma X4 of |23j); 
we also refer to [23] for a general discussion. The following two lemmas (see Chapter 5 
of [25]) have been used by many authors (for instance, see Lemmas 1, 2 and 3 of [26] or 
Lemmas 2.3 and 2.4 of [23) to control the nonlinear terms. 

Lemma 3.3 Let s > 0, / G C [s]+1 (M) with f(0) = 0. Then for any u G H s (R n ) n 
L°°(IR n ) ; we have f(u) G # s (IR n ) n L°°(R n ). Moreover there is some constant A(M) 
depending on M (and s) such that for all u e H s (R n ) n L°°(M n ) with \\u\\ L o° < M 

\\f(u)\\ s <A(M)\\u\\ s . 

Lemma 3.4 Let s > 0, f G C^ +1 (M). Then for any M > there is some constant 
B(M) depending on M (and s) such that for all m,u 2 G H s (R n ) n L°°(R n ) with 
||wi||l°° < M, IIM2IU 00 < ond || ||s < M, || ^2 1| s < M we have 

\\f{u x )-f{u 2 )\\ s <B{M)\\u x -u 2 \\ s . 

In our case the nonlinearities are of the form 
BF BF 

^-(|V«;| 2 ) = 2w x F'(\Vw\ 2 ), |^(|V^| 2 ) = 2w y F'(\Vw\ 2 ) 

where F' denotes the derivative of F. It follows from repeated applications of Lemma 
13.21 that for the above terms Lemmas 13.31 and 13.41 take the following forms: 

Lemma 3.5 Let s > 1, F G C [sl+1 (K)- Then for any w G X s , we have 

BF BF 

t^(\Vw\ 2 ) G H CR 2 ), ^±-(|Vw| 2 ) G H 8 -\R 2 ). 

BW X OWy 

Moreover there is some constant A(M) depending on M (and s) such that for allw G X s 
with || |Vw| 1 1 1,00 < M 
BF 

" ' ( I | 2 ) || s-l < j4.(M)||u>|| s , 



Bw x 
BF 

\\^(\Vw\ 2 )\\ s ^ < A(M)\\w\ 
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Lemma 3.6 Let s > 1, F G Then for any M > there is some constant 

B(M) depending on M (and s) such that for all w 1 ,W2 G X s with — M, 

\\w2\\x a < M we have 

BF BF 
\\— (\V Wl \ 2 ) - — (|V«; 2 | 2 )IU < B(M)\\ Wl - w 2 \\„ 

BF BF 
H^dV^I 2 ) - — (|V W2 | 2 )|U < B{M)\\ Wl -w 2 \\ s . 

We want to emphasize that, while the constants B(M) in Lemmas I3.4l and 13.61 are indeed 
local Lipschitz constants depending on the X s -norms, the constants A(M) in Lemmas 
13.31 and 13.51 depend only on the L°°-norms. This property of A(M) will be used in the 
proof of Lemma 13.91 which characterizes the type of blow-up. 

When s > 2 we have the following local well posedness result. 

Theorem 3.7 Suppose s > 2 and the decay rate r in M.3\) satisfies r > 2. For 
tp,ip G H S (R 2 ), there is some T > such that the Cauchy problem M.l\) - U~^ is well 
posed with solution w = w(x,y,t) in C 2 ([0, T], H S (R 2 )) . 

Proof. Let w G H S (R 2 ). For s > 2, by the Sobolev Embedding Theorem we have 
|Vw| G L°°(R 2 ). Thus X s = H S (M. 2 ) and the norm ||u>||x a can be replaced by the 
equivalent H s norm ||u>|| s . Since f ll.ip - fll.2p is formally equivalent to f l3.ip -f l3T2|) . we will 
use the well-known existence-uniqueness (Picard-Lindelof) theorem for Banach space 
valued systems of ordinary differential equations (for instance, see Theorem 5.1.1 of 
[28J). Obviously, all we need is to show that the operator K of (13. 3 p is locally Lipschitz 
on X s . We first show that K maps X s into X s . We estimate the convolution as 

\\(3*u\\ s = n(i + \er /2 mm\\ < cm + m (s - r)/2 um = ^imu, 

where we have used inequality (11. 3p . By Lemma [3.51 for || |Vu>| ||i°o < M 

BF \ M lln BF „ BF „ 



Bw x J x Bw x Bw x 

< CA(M) \\w\\ s+2 - r < CA(M)\\w\\ a 

where we have used r > 2. The same holds for the term ( f3 * 4^- J and 

\\Kw\\ s < CA(M)\\w\\ s+2 _ r < CA(M)\\w\\ s . (3.5) 
Similarly, for Wi,w 2 G X s with ||iOi|| s < M and \\w 2 \\ s < M, by Lemma [3.61 

< CB(M)\\wi-w 2 \\ s . 

As above, the same holds for the term f (3 * 4^- ) . So, K is locally Lipschitz on X s and 

thus the local well posedness of the Cauchy problem is established. ■ 

When r > 3 in fll.3p . the extra regularizing effect of /3 allows us to improve the 
result in Theorem 13.71 to the case of s > 1 . 
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Theorem 3.8 Suppose s > 1 and the decay rate r in U.3\) satisfies r > 3. For 
(p,i/j £ X s , there is some T > sttc/i t/iat the Cauchy problem U. ip -z TOj) is u>eW posed 
with solution w = w(x,y,t) in C 2 ([0,T],X s ) . 

Proof. Similar to the proof of Theorem 13.71 it suffices to show that the map K given 
in (I3.3p is locally Lipschitz on X s . Recall that ||-frw7||x s = ||-Kw|| s + + 
|| (Kw) y \\Loo. The term ||iCiu|| s can be estimated by ||io|| s as above. For e = r — 3 > 
we have 

\\(Kw) x \\ L oo < C\\(Kw) x \\ 1+e < C\\Kw\\ 2+e < C\\Kw\\ s+1+c 

< CA(M)\\w\\ s+3+e - r = CA(M)\\w\\ s . (3.6) 

where we have used (13. 5p and the Sobolev Embedding Theorem. The same holds for 
(Kw) y and a similar estimate as in the proof of Theorem 13.71 shows that K is locally 
Lipschitz on X s . ■ 

The solution of ( ll.ip - (jl.2p can be extended to a maximal interval [0,T max ) where 
finite T max is characterized by the blow up condition 

limsup (\\w(t)\\ X s + \\w t (t)\\x') = oo. 

n x 

Obviously T max = oo, i.e. there is a global solution, if and only if for any T < oo 
limsup (||w(t)||x= + \\wt(t)\\ X s) < oo. 

The lemma below characterizes the type of blow-up; namely blow-up occurs in the 
L°°-norm of |Vw|. 

Lemma 3.9 Suppose that the conditions of Theorem \3. 7| or Theorem \3. 8\ hold. Then 
T max = oo ; i.e. there is a global solution of the Cauchy problem M.l\) - [1~^) . if and only 
if for any T > 

limsup( 1 1^(^)1 1 L°o + \\w y (t)\\ L ™) < oo. 

t->T- 

Proof. Since 

II^WIU 00 + ll w y(*)IU°° < IkWIIx-, 
it suffices to prove that if the solution exists for t £ [0, T) and 11^(^)1^°° + ||wy(£)|| l°° < 
M for all < t < T then both ||w(t)||x s and ||wi(t)||x s stay bounded. Integrating 
equation (II. ip twice and calculating the resulting double integral as an iterated integral, 
we obtain 

w(t) =<^ + t^+ / it - t){Kw){t)cLt, (3.7) 



J o 

w t (t)=i/>+ / (Kw)(r)dT. (3.8) 
Jo 

But, by (E3D, \\(Kw)(t)\\ s < CA(M)\\w(r)\\ s+ 2-r < CA(M)\\w(t)\\ s where the constant 
A(M) depends only on M. Hence 

\Ht)\\s + \\wt(t)\\. < Ms + (1 + T)||^|| s + (1 + T)CA(M) [ \\w(r)\\ s dr, 
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and Gronwall's Lemma gives 

IK«)II« + \H(t)\U < (Ms + (1 + T)U\\ s )e^ CA ^ T (3.9) 

for all t G [0,T). We now estimate ||u>fcc(£)||.L«>. The estimate for ||i%(£)||l°° follows 
similarly. In the case of Theorem 13.71 (where s > 2), by the Sobolev Embedding 
Theorem, 



||wta(t)IU°° < C\\w t ( 
so that (13.91) applies. In the case of Theorem 13.81 (where r > 3), from (I3.8P 

\\w tx (t)\\ L oo < \\^ x \\ L oc + / \\(Kw) x (t)\\ l °o dr. (3.10) 



o 



By m> 

|| {Kw) x (r) |Uco <CA{M) II w{t)\\ 3 , 

and by (13~91 

lh(r)|| s <(||^|| s + (l+T)||^|| s )e( 1+ ^^. 
Finally, plugging into (13.101) we obtain the required estimate 

\\w tx (t)\\L°° < UAloo + TCA{M)(\\<p\\. + (1 + T)|H| s ) e ( 1+T ^ M ) T , 
which corresponds to the case s > 1, r > 3. ■ 

4. Conservation of Energy and Global Existence 

In the present section we will prove that locally well defined solutions can be extended 
to the entire time. 

In the study of global existence of solutions the conservation of energy plays a 
key role. First, time invariance of the energy functional will be shown. To this end, 
we define an unbounded linear operator R on H S {R 2 ) as Ru = J 7-1 ((/3(£))~^(0) 

where J 7-1 denotes the inverse Fourier transform and /3(£) is defined in (II. 3p . Obviously, 
RPu = J 7-1 ^(/3(0)~^(£)) f° r a rea l number p with u in the domain of RP. On the 
other hand, R~ 2 u = j3 * u and formally (II. ip can be rewritten as 

*v = m + m . <4,) 



dw x J x \dw y , v 

Here we have used the fact that convolution commutes with differentiation in the 
distribution sense, i.e. ((3 * u) x = (3 * u x . 

Lemma 4.1 Suppose that the conditions of Theorem \3.7\ or Theorem \3.8\ hold and the 
solution of the Cauchy problem /TO)- /TO) exists in C 2 ([0, T) , X s ). If Rip E L 2 (R 2 ), 
then Rw t {t) e L 2 (R 2 ) for allt e [0,T). Moreover, ifR<p e L 2 (R 2 ), then Rw(t) G L 2 (M 2 ) 
for allt G [0,T). 
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Proof. Formally, from ( 13. 8 p we have 

Rwtit) = Rip + [ {RKw)(r)dT . (4.2) 
Jo 



Note that 



/ dF \ ( dF 
RKw = [a * — — + \ a* 



® w x J x V 9w y 

where a(£) = (/3(^)) 1 / 2 . Then similar to the derivation of (I3.5P (replacing j3 by a and 
hence r by r/2) we get 

||(it;^)(T)|| s+f _ 2 <c|| w (T)|| s . 

Since either (s > 2 and r > 2) or (s > 1 and r > 3), in both cases we have s + 1 — 2 > 0. 
Thus the right-hand side of ( 14.21) belongs to L 2 (IR 2 ) and the conclusion follows. The 
second statement follows similarly from (13. 7p . ■ 

Lemma 4.2 Suppose that the solution of the Cauchy problem ( fi. ij) -( TOj) exists on some 
interval [0,T). If Rip G L 2 (R 2 ) and the function F(\V(p\ 2 ) belongs to L 1 ^ 2 ), then for 
any t G [0, T) the energy 

E(t) = -\\Rw t (t)\\ 2 + I F(\Vw(t)\ 2 )dxdy (4.3) 

is constant in [0, T) . 

Proof. By Lemma \A. 11 Rwt(t) G L 2 (IR 2 ). Multiplying (14.11) by Wt, integrating in x and 
y, and using Parseval's identity we get 

° = It L l^' 1 ^^ + J R2 F\\Vw{t)\ 2 )^{\Vw{t)\ 2 )dxdy 

= JtL Q (iM(t))2 + F ^ Vw ^) dxdy 
which implies the conservation of energy. ■ 

The main result of this section is the following theorem. 

Theorem 4.3 Suppose s > 1 and the decay rate r in U.3\) satisfies r > 4. Let 
(p,ip G X s , Rip G L 2 {R 2 ) and F(|V<^| 2 ) G L 1 ^ 2 ). If there is some k > so that 
F(u) > —ku for all u > 0, then the Cauchy problem / ti.ip -/ TOp has a global solution in 
C 2 ([0,oo),X*). 

Proof. By Theorem 13.81 the Cauchy problem is locally well-posed. Assume w G 
C 2 ([0,T),X s ) for some T > 0. Since F(u) > -ku, for all t G [0,T) we have 

||iM(t)|| 2 = 2E(0) -2 [ F(\Vw(t)\ 2 )dxdy, 

<2E{0) + 2k / \Vw{t)\ 2 dxdy, (4.4) 

./R 2 
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where E(0) is the initial energy. On the other hand we have 

HiMt)ii 2 = / m)r l \wt(t,t)\ 2 d£, 

Jr 2 

>C-' jjl + \^\wt(U)\ 2 ^, 

= C^ml, (4.5) 
where (jl.3p is used. Combining (I4.4p and (14. 5 p 

|lK0lll<2|Kt)|| s IK(t)|| s 

< lKt)||| + IK(t)||| 

< \Ht)\\\+C\\RvH(t)\\ 2 

< \\w{t)\\ 2 r + 2C{E{0) + k\\\Vw{t)\\\ 2 ) 

< 2CE{0) + (1 + 2CA;)|K*)|||, 

where |||Vw(t)||| < ||w(t)||i < ||w(t)||r is used. Gronwall's lemma implies that ||iu(t)|| r 
stays bounded in [0, T). As r > 4 we have | — 1 > 1 and the Sobolev Embedding 
Theorem implies 

lllv^Hi^^iiiv^iiir-x^iKOIIi- 

We conclude that |||V^(t)|(t)|| LOO also stays bounded in [0,T). By Lemma 1331 this 
implies a global solution. ■ 

5. Blow up 

In this section a blow-up result for (ll.ip -( TL~2]) will be presented. The following lemma, 
based on the idea of Levine , will be used to prove blow up of solutions in finite time 
for certain nonlinearities and initial data. 

Lemma 5.1 Suppose that T-L(t), t > 0, is a positive, twice differentiable function 
satisfying U"U - (1 + u){W) 2 > where v > 0. If H{0) > and H'(0) > 0, then 
Hit) -> oo ast->ti for some h < U(0) / 'vW '(0) . 

Theorem 5.2 Suppose that the solution, w, of the Cauchy problem UJ]) - l[Tlfy exists, 
Rip, Rip e L 2 (R 2 ) and F(|Vy?| 2 ) e L^IR 2 ). If there exists a positive number v such 
that 

uF\u) <(l + 2v)F{u) forallu>0, (5.1) 

and 

£(0) = i||i#|| 2 + / F(\V<p\ 2 )dxdy<0, 
z Jr 2 

then the solution, w, of the Cauchy problem U. ip -z TOj) blows up in finite time. 
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Proof. We assume that the global solution to ( ll.ip - (ll.2p exists. Then, by Lemma 14.11 
Rw(t), Rw t {t) e L 2 (R 2 ) for all t > 0. Let H{t) = \\Rw(t)\\ 2 + b(t + t ) 2 where b and t 
are positive constants to be determined later. Then we have 

W = 2(Rw t , Rw) + 2b(t + t ) 

U" = 2\\Rw t \\ 2 + 2(Rw tt , Rw) + 2b. 

Note that %'(0) = 2(R<p,Rip) + 2bt > for sufficiently large t . Using the inequality 
uF'(u) < (1 + 2v)F(u) together with (IQl) and (14T51) we have 

(Rw tt ,Rw) = (R 2 w t t,w) 

= -2 1 \Vw\ 2 F'{\Vw\ 2 )dxdy 
Jr 2 

> -2(l + 2v) / F(\Vw\ 2 )dxdy 
Jr 2 

= (l + 2u)(\\Rw t \\ 2 -2E(0)), 

so that 

H" > 4(1 + v)\\Rw t \\ 2 -4(1 + 2^)^(0) + 2b. 
On the other hand, using 2ab < a 2 + b 2 and Cauchy-Schwarz inequalities we have 
(Uf =4[(Rw,Rw t )+b(t + t )] 2 

<4(||i?w|| 2 \\Rw t \\ 2 + 2b(t + t )\\Rw\\ \\Rw t \\+b 2 (t + t ) 2 ) 

< 4 (\\Rw\\ 2 \\Rw t \\ 2 + b\\Rw\\ 2 + b\\Rw t \\ 2 (t + t ) 2 + b 2 (t + t ) 2 ) . 

Thus 

H"H- {l + v)(H') 2 

> (4(1 + ^)||i?^|| 2 - 4(1 + 2v)E{0) + 2b) (\\Rw\\ 2 + b(t + t ) 2 ) 

- 4(1 + v) (\\Rw\\ 2 \\Rw t \\ 2 + b\\Rw\\ 2 + foll^H 2 ^ + t ) 2 + b 2 (t + t ) 2 ) 

= -2(l + 2u)(b + 2E(0))H. 

Now if we choose b < —2E(0), this gives 

U"(t)H(t) - (1 + v) {H\t)f > 0. 

According to the Blow-up Lemma \5.1\ this implies that R{t), and thus ||i?w(t)|| 2 blows 
up in finite time. ■ 

Consider a typical nonlinearity of the form F(u) = au q with q > 0. When a > 0, 
Theorem 14.31 will apply and a global solution exists (for suitable s, r and initial data). 
When a < 0, the blow-up condition (15. ip of Theorem 15.21 holds if and only if q > 1. 
This observation says that the global existence result of Theorem 14.31 is essentially sharp 
compared to the blow-up result. 

Finally, we conclude with a short discussion on the condition E(0) < 0. If F(u) > 
for all u > 0, then by Theorem 14. 3[ there is a global solution. If F(u) is negative on some 
interval J, we can choose ip with support in / so that f R2 F(\V(p\ 2 )dxdy < 0. Hence, 
when ip = or Rip is sufficiently small we get E(0) < 0. This also shows that blow up 
may occur even for small initial data. 
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6. The Anisotropic Case 

In this section we will consider the more general nonlinear term of the form F(w x , w y ) 
with F(0,0) = rather than the isotropic form F(|Vw| 2 ). Such a form appears as the 
strain energy function of anisotropic materials [18]. With minor modifications on the 
assumptions, all the results of Sections 2-5 can be generalized. Following the layout 
of the manuscript we will pinpoint out these modifications and briefly explain how the 
proofs will change accordingly if we replace F(\ Vw\ 2 ) by F(w x , w y ) = F(Vw) where the 
symbol ~ is employed to distinguish the anisotropic form of the strain energy function. 



Local Existence and Uniqueness of Solutions 

The main step in Section 2 is to show that the map K of (13.31) is locally Lipschitz on 
the Banach space X s . To deal with the more general nonlinearity F(Vw) we will use 
the vector versions of Lemmas 13.31 and 13.41 |25j. For a vector function U = (ui,u 2 ) the 
notation \\U\\ = \\ux\\ + ||w 2 || will be employed where |.| denotes a given norm. 

Lemma 6.1 Let s > 0, h G C W+1 (K 2 ) with h(0) = 0. Then for any U = (ui,u 2 ) G 
(iP(R 2 )nL°°(IR 2 )) 2 ; we have h(U) G F s (M 2 )nL°°(lR 2 ). Moreover there is some constant 
A(M) depending on M such that for all U G (H S (M. 2 ) n L°°(R 2 )) 2 with \\U\\ LOO < M 

\\h(U)\\ s <A(M)\\U\\ s . 

Lemma 6.2 Let s > 0, h G C W+1 (K 2 ). Then for any M > there is some constant 
B(M) such that for all U U U 2 G (H s (R 2 )nL°°(R 2 )) 2 with ||C/i|| LOO < M, \\U 2 \\ Loa < M 
and \\Ui\\ s <M, \\U 2 \\ S < M we have 

\\h{U 1 )-h{U 2 )\\ s <B{M)\\U 1 -U 2 \\ s . 

Through Lemmas 16. II and 16.21 we obtain the estimates of Lemmas 13. 51 and l 3.6l for F(Vw). 
With these new estimates, the proofs of Theorems 13.71 and 13.81 on local well posedness 
and of Lemma 13.91 about the blow up criterion follow exactly the same way for the 
anisotropic form. 

Conservation of Energy and Global Existence 

Due to the new estimates on the map K of (13.31) . Lemma \A . 1 1 holds also for the general 
case. In Lemma 14.21 the energy must be replaced by 

E{t) = h\Rw t (t)\\ 2 + / F(Vw(t))dxdy. 

The proof of the energy identity is the same with the obvious modification 




r r I oF oF \ 

I F(Vw(t))dxdy = / - — w xt + - — w yt dxdy 
Jr2 J R 2 \ dw x dw y ] 

w t dxdy. 



y 

'dF 



dlUy 
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In Theorem I4.3[ we change the assumption on the nonlinearity as F(U) > — k\U\ 2 for 
all U € M 2 . Then the estimate f !4.4j) still holds for the anisotropic case and thus the rest 
of the proof follows. We note that our new assumption F(U) > — k\U\ 2 reduces to the 
previous condition F(u) > —kuoi Theorem 14.31 when the nonlinearity is of the form 
F = F(\Vw\ 2 ). 

Blow up 

The blow-up condition 

uF\u) < (1 + 2v)F{u) for all u > 0, 
of Theorem 15.21 takes the form 

U ■ VF(U) < 2(1 + 2v)F(U) for all !7g1 2 

for the general nonlinearity F(Vw). With this new blow-up condition all the steps in 
the proof of Theorem 15.21 are still valid for the general case except for the modification 
below in the estimate for the term (Rw u , Rw): 

(Rwu, Rw) = (R 2 w tt , w) 

Vw • VF(Vw)dxdy 

> -2(l+2i/) / F(Vw)dxdy 

</R2 

= (l + 2u)(\\Rw t \\ 2 -2E(0)). 

We note that the new blow-up condition reduces to the blow-up condition of Theorem 
15.21 for the isotropic form F(|Vw| 2 ). 
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